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I. INTRODUCTION 

Spectacular advances in bio- and nanotechnology make 
it possible, not only to measure or observe, but also to 
manipulate and construct objects at a very small scale. 
At the same time there is growing interest in techniques 
which can add functionality. In particular, the devel- 
opment of molecular engines is a theme which has re- 
ceived great attention over the last two decades. The ap- 
pearance of fluctuations in small systems has led to new 
concepts for characterizing or operating such devices, ex- 
ploiting rather than fighting these very same fluctuations. 

These so-called Brownian motors [TJ [5] have an ad- 
ditional theoretical interest through their relation with 
the old issue of Maxwell demons and the second law of 
thermodynamics. In return, this theoretical connection 
allows to make statements on the efficiency of such en- 
gines [3] or to transform them from engines into mini- 
refrigerators [H [5] . Most of the studies on Brownian mo- 
tors start with an ad hoc separation of systematic and 
noise terms, based on linear Langevin equations. This 
approach however offers little insight into the origin of 
the rectification of random fluctuations. As pointed out 
by van Kampen [B], the rectification of nonlinear fluc- 
tuations cannot be addressed starting from the standard 
Langevin description with additive Gaussian white noise. 
In [71 [5] a theoretical and numerical study of a thermal 
engine is presented in which rectification arises at the 
level of nonlinear response. The analysis therein starts 
from a microscopic description based on Newton's laws of 
motion. There is another related distinct feature of the 
model: the asymmetry of the thermal engine lies in the 
geometry of the motor itself, in contrast to the asymme- 
try imposed by the application of an external potential, 
appearing in the so-called flashing and rocking ratchet 
models. 

The characteristic properties of the engine, such as 
the friction coefficient, the speed and the diffusion co- 
efficient, are calculated exactly in [71 IH] and are found to 
be in excellent agreement with the results from hard disk 
molecular dynamics. However, the results are reported 
in dimcnsionless units, in part due to the fact that the 
analysis was, for reasons of simplicity and for compari- 
son with molecular dynamics, limited to the case of two 
dimensions. In view of the technological interest of mo- 
tors in bio- and nanotechnology, we report here a full and 
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FIG. 1: The two-reservoir model of the thermal engine: two 
solid objects are confined in separate containers (scaled down 
for illustration purposes) that contain gases at temperatures 
Ti and Tb. The objects are assembled with a rigid connection. 
The ensemble can move freely along the z-axis. 

detailed analysis of the three-dimensional version. 

This paper is organized as follows. First, the model, 
notations, and working hypothesis are introduced in sec- 
tion |nj The calculation method, based on the kinetic 
theory of gases, is presented in section |III| with a dis- 
cussion of the analytical solution following in section |IV| 
Finally, in section [V] we report and discuss the results for 
the friction coefficient, diffusion coefficient, and drift ve- 
locity as a function of shape and present estimates based 
on physically realistic parameter values. 

II. THE MODEL 

The model presented in [7J reproduces in a simplified 
way the principle ingredients of Feynman's ratchet and 
pawl mechanism [9]: a temperature difference between 
two reservoirs and the presence in at least one reservoir 
of an asymmetric object. The construction, extended to 
the case of three spatial dimensions, is as follows. We 
consider any number of reservoirs (denoted by index i), 
each containing a gas at equilibrium at a temperature Ti. 
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Fig. [T] gives a schematic picture of the two-reservoir sys- 
tem. Solid objects with no internal degrees of freedom, 
called 'motor units', are located inside the containers. 
These objects are coupled rigidly to each other, so that 
the motor moves as a single entity, with total mass M, 
along a given straight axis, corresponding to its single 
translational degree of freedom. For simplicity, we dis- 
regard any rotational degree of freedom (for a detailed 
discussion of this case, see [5]). 

Due to collisions with the gas particles (mass to), the 
motor will change its velocity V(t) in the course of time t. 
The statistics of these collisions can be described under 
the assumption of molecular chaos, which is valid when 
the gases are in the high Knudscn number regime and 
the containers are large enough to avoid acoustic and 
other boundary effects. In addition, the shape of the 
units' surfaces must be such that no re-collisions with 
the motor occur, namely for convex and closed shapes. 
With these assumptions, the precollisional velocities are 
random and uncorrelated. Hence the time evolution of 
the probability P(V, t) that the motor has speed V at 
time t can be described by a master equation: 



dP(V,t) 
di 



= / dV [W{V\V')P{V',t) - W(V'\V)P(V,t)] . 

(1) 

Here W(V|V) represents the transition probability per 
unit time for the motor to change its speed from V' to 
V. 



III. KINETIC THEORY 



We will also suppose that the collision is described in 
terms of a (short-range) central force, implying that the 
component of the momentum of the gas particle along 
any direction tangential to the surface of the motor is 
conserved. The orientation of the tangent plane to the 




FIG. 2: (Color online) The orientation of an infinitesimal sur- 
face element dS is represented by an outer-pointing unit nor- 
mal vector e± and determined by the spherical coordinates ip 
and 6. The polar angle 8 is measured from the z-axis, which 
is chosen as the free direction of movement of the motor, x, y 
complete the Cartesian coordinate system and the azimuthal 
angle ip starts from the x-axis. Two orthogonal unit vectors 
ei jM and e*2,n determine the plane that is tangent to the motor 
unit in dS. 

motor surface is determined uniquely by a normal out- 
ward vector on an infinitesimal element dS of the surface 
at the point of collision. In spherical coordinates, this 
normal vector is given by 



In this section we study the collisions of gas particles 
from either temperature reservoir with a motor part and 
derive the resulting total transition probability IF(V|V') 
for the motor to change speed from V' to V. We intro- 
duce a Cartesian coordinate system (a;, y, z) where the 
z-axis points along the free direction of movement of the 
motor. 



A. Conservation rules 

A gas particle will, upon collision with a motor unit, 
undergo an instantaneous change of velocity from v' — 
(v' x ,v' y ,v z ) before collision to v — {v x ,v y ,v z ) afterwards. 
Due to conservation of momentum along the free z- 
direction, one has: 



/ + MV' = rnv z + MV. 



(2) 



In addition, when the collision is perfectly elastic, the 
total energy is conserved: 



-MV' 2 
2 



1 ,2 

1 



MV 



1 ,2 

2 mv v 

1 



1 ,2 

2 z 



-mil, 



1 2 

2 mv v 



-mu, 



(3) 



e_L I = (sin 9 cos (p, sin 9 sin ip, cos ( 



(4) 



with 9 the polar angle from the z-axis (0 ^ 9 ^ 7r) and 
tp the azimuthal angle in the £j/-plane from the x-axis 
(0 < ip < 2n), cf. Fig. [2} We also introduce two mutually 
perpendicular unit vectors within the tangent plane: 

ei, M = (cos # cost/?, cos # sin — sinfl), (5) 
6*2,11 = (-sin (p, cos (f, 0), (6) 

so that we can write the conservation of tangential mo- 
mentum as 



V ■ ei jM = v ■ ei,„, 
v' ■ e 2 ,n = v ■ e 2 .,. 



(7) 
(8) 



Solving the conservation rules [Eqs. ( 2|3|7|8 )] for 

V, v x ,v yi v z leads to the following expression for the post- 
collisional speed V of the motor in terms of the precolli- 
sional speeds: 



V = V 



2^cos : 



?9 



1 + g cos 2 9 



I (v' x tan 9 cos ip - 



(9) 

v' y tan 9 sin ip + v' z — V 1 ). 



(10) 
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B. Transition probability 

The motor is subject to random collisions by gas par- 
ticles form the different reservoirs i. The particle den- 
sity and velocity distribution in reservoir i are denoted 
by pi and <fri(v x ,v y ,v z ), respectively. The contribution 
dWi(V\V) to the total transition probability W(V\V), 
coming from collisions on an infinitesimal surface clement 
dSi of the motor unit in reservoir i, can then be found by 
considering the number of gas particles that collide with 
this surface element dSi in a unit time step: 



dWi(V\V) = dSi 



+00 



+ OC 



dv' y 



dv' 



x H[(V> - v>) ■ e ± ]\(V' - v') ■ e ± | Pi <M<,<,<) 
x S[V-V 

- B(0) (tan cos ip v' x + tan sin ip v' y +v' z -V')]. 

(11) 

Here H represents the Heaviside function, S the Dirac 
distribution, and 



om ros 2 n 

B{6)= M 



1 



M 



cos 2 6 



(12) 



The Dirac delta distribution selects those particles that 
produce the required postcollisional speed V , following 
the collision rules [Eq. ( 10 1] . Assuming that the veloc- 
ity distribution of the gas particles is Maxwellian at the 
reservoir temperature Tj, 



4>i(V X ,Vy,V Z ) 



2Trk B Ti 



3/2 



exp 



m{vl 



2k B Ti 



(13) 

the integrals over v' x , v' , v' z in Eq. (Ill can be calculated 
explicitly. The total transition probability for the motor 
to change velocity from V' to V in a unit time, is then 
found by integrating dWj(F|V) over the surface Si of 
each motor part and summing over all the reservoirs i: 



W{V\V') = W> 



2-Kk B Ti 



x ^(V - V')H[V - V'] J 

■(V'-V)H[V'-V] f 
Js 



dSi 



cos e>o 



AI 



/Si,cos6»<0 



mcosf 



cos 9 



xexp 



III 2-t--> 1- M 



2k B Ti 



cos 2 0[ ^' + -(1 + 



mcos^ 



-XV -V) 



with P eq (V) the Maxwell Boltzmann distribution for the 
speed V of the motor unit (at the temperature of the 
reservoir(s)). This is in agreement with the general prin- 
ciple of detailed balance in a system at equilibrium |10j . 



IV. SOLUTION METHOD 

We follow the method of [S] to solve the master equa- 
tion [Eq. 0] for the moments of the motor velocity, 



(V n ) = / P(V,t)V n dV. 



(16) 



This method is based on the van Kampen 1 /Q expansion 



[6J- 



A. Solution of the master equation 

It is convenient to scale the motor velocity V to a di- 
mensionless variable 



X 



M 



k B T e g 



-V, 



(17) 



with the effective temperature T c g to be determined self- 
consistently from the condition (X 2 ) = 1 in steady state 
operation. We can expand the integrand in Eq. ([!]) in a 
Taylor series about X': 



v (-1)" d" 



Here the 'jump moments' are given by 



J n {X)= J AX n W(X;AX)dAX, 



(19) 



with W(X'; AX) = W(X\X') and AX = X'-X. Using 
Eq. ( 18 1 a coupled set of equations for the time evolution 



of the moments (X n ) can then be constructed: 



d(X n ) 
dt 



k=l 



(20) 



with (?) the binomial coefficients. 



We note that in the case of a single reservoir or for mul- 
tiple reservoirs at the same temperature, the transition 
probability satisfies the following relation: 

W{V\V')P eq {V) = W{-V'\ - V)P eq (-V), (15) 



The exact expression for the jump moments J n (X) is 
obtained by integration over AX in Eq. (19). In terms 
( 14 bf parabolic cylinder functions [IT] . 

D„(z) = (exp[-z 2 /4]/r[-n]) 



m / exp[-zx- x 2 /2]x~ n ~ 1 dx (for n < 0), (21) 
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and the Gamma function (T), the result reads: 

MX) 



2™ r JIM 
V2n \ V m 




< I dSi I (■(>> f 



M 1 



m cos ( 



x exp 
x D_„_ 2 



AM Ti 
1 /m / T rS 



2 V M V T, 



A cos (9 



(22) 
(23) 

(24) 
(25) 
(26) 



An exact solution of Eq. (20 1 is not available. We 



therefore turn to a perturbational approach in terms of 
the parameter 



= ^Jm/M. 



(27) 



This is consistent with the observation that the mass M 
of the motor is expected to be much larger than the mass 
m of the gas particles. Even for a motor with dimensions 
of nanometers operating in a gaseous environment, e is 
of order of 10 -3 . 

The expansion of the parabolic cylinder functions is 
given by 

2-/ 2 I> + 2] D_„_ 2 (z) = - V2T[^]z 

2n + 3 n + 2 n 2n + 3 T , r n + 3 1 3 



4 3 
4n 2 + 12n+ 11 
96 



6\/2 ~ L 2 
r[(2 + n)/2]z 4 + 0(z 5 ). (28) 



Introducing a scaled time r = e 2 t, we find the following 
expansion for the equation of the first moment: 



!f 2 / 2 {x 

J fi ff V 7T 




+ 0(e 4 ). (29) 



The geometry of the motor is contained in the shape 
factors cr nji , defined as: 



a n ,i = / dSi cos r ' 
I Si 



(30) 



At this point we remark that the azimuthal angle ip has 
dropped out, and that the geometric dependency is de- 
termined only by the (polar) angle 6 between the surface 
and the direction of movement. Also note that the term 



me 1 in Eq. ( |29| is zero by application of Gauss' theo- 
rem: 



G 1, 



dSi cos 8 — / dSie± ■ e z 
Si Js, 



dVi(V-e z ) = 0, 
(31) 

where the latter integral is over the interior volume Vi 
of a motor part. This is consistent with the fact that 
there is no net macroscopic force acting on the motor. 
The net motion that will be revealed below is the effect 
of fluctuations only. 

Similarly, for the equation of the second moment, one 
finds the following expansion: 



d(X 2 




-2e' 



Jeff o J-i 



Linear relaxation 



+0(e 3 ). 



(32) 



To order e°, Eq. (29 1 reduces to a linear relaxation 



law Mdt(V) = — 7(V)with 7 = ^ 7, the sum of linear 
friction coefficients 7$ of each part of the object: 



4p, 



k B Tim 



2tt 



-02,; 



(33) 



where Qi = mpi is the mass density of the gas, = 
y/8kBTi/ (nm) is the mean gas velocity, and 02, i is a 
geometric factor. 



C. Nonlinearity: Steady state directed motion 

If a constant temperature difference between the reser- 
voirs can be maintained for a time longer than the relax- 
ation time M/7, the probability distribution will relax 
to its steady state value. Restricting ourselves to the 
first two moments, we turn to the steady state solution 



of Eqs. (29 1 and (32 1. First, from Eq. (32) we determine 



the effective temperature T c g, which was defined earlier 
by the condition (X 2 ) = 1. To lowest order, e°, we find 
that T c g is the weighted average of the reservoir temper- 
atures: 



Teg — 



Ei7i 



i Pi°%iT\ 



(34) 
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TABLE I: The geometric moments 02 and 03 for some basic three-dimensional shapes. The lowest orders 02 and 0-3 (01 = 0) 
are tabulated in terms of the total surface area S of the geometry. For illustrations of the shapes and their parameters, see 

Fig-j 

Shape &2/S o-s/S Surface S 

Disk 1 2nr 2 

Blade 1 2lw 

Sphere 1/3 4?rr 2 

Cone sina sina(sina — 1) 7rr 2 (l + esc a) 

Pyramid sina sin a(sin a — 1) \nr 2 cot - (1 + esc a) 

Spherical cap cos 2 °+ 5 c l° °+ 6 — = nr 2 (3 + cos a) /{I + cos a) 

1 3cosa+9 cos2a+4cosa-5 \ 1 // \ ' / 

Spherical cone 3 sin 3 Q -2cos 3 o,+ 2 2 sin 4 a +cos 4 a -i yrr 2 esc a(l + 2 tan(a/2)) 

r Jsin»-Dcosa+D 4 cos a — A sin a — 4 v \ f // 




FIG. 3: Examples of simple three-dimensional geometries that can be treated analytically. The parameters that determine the 
relative areas of the surfaces are indicated in the illustrations. 



Solving Eq. ( 29 1 to order e° leads to a zero average drift 
speed (X) = 0. Rectification of the thermal fluctuations 
occurs at higher levels of the expansion and it is necessary 
to include nonlinear terms. Solving Eq. (29 1 up to order 



e gives us an expression to lowest order for the average 
drift speed of the motor: 



(V) 



in lirk B T eS J2iPi(^J. 



8M 



-In; 12. 12 j PiPjV3,i<72,j(Ti - Tj)y/Tj 



M V 8M 



(12iPiVTi(T2,: 



V. RESULTS AND DISCUSSION 



Friction and diffusion coefficients 



(35) 




10 20 30 40 50 60 70 80 90 

Opening angle a 

FIG. 4: Diffusion coefficients of objects with different geom- 
etry but identical cross section irr 2 (r — 100 nm) in highly 
diluted argon gas (10 19 m~ 3 particles, mean free path 2.3 /zm) 
at 299 K temperature. The geometries are (a) a disk, (b) a 
sphere, (c) a cone, (d) a spherical cap, and (e) a spherical 
cone. See Fig. [3] for illustrations. For (c,d,e) the shapes and 
hence the diffusion coefficients depend on an opening angle a. 



Although not directly related to the main topic of this friction coefficient of corresponding basic shapes. The 
work, we briefly pause to discuss the new result for the result for a sp h e rical shape with radius r, 
linear friction coefficient, given in Eq. (33 1. Together 

with results from Table |T] for the geometric factors er 2 4 2 

this result provides the explicit expression for the linear 7 — ~^ 7Tr Qi v i> 

(36) 
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is in agreement with the result found in [12 . In combi- 
nation with the Einstein relation D = ksT/j, we also 
obtain the explicit formulas for the corresponding diffu- 
sion coefficients D. As an illustration, numerical values 
are given in Fig. [I] for various shapes of cross section nr , 
r = 100 nm. To be concrete, we consider highly diluted 
argon gas (density 10 19 m -3 ) leading to diffusion coef- 
ficients of the order of 3 x 10~ 4 m 2 /s, and correspond- 
ing friction coefficients of order 1.5 x 10~ 17 Ns/m. As 
expected, the conical shapes have higher diffusion coeffi- 
cients and lower friction as one considers smaller opening 
angles a. 



B. Equilibrium 

When the reservoirs are all at the same temperature, 
Ti = T, we find T c g — T, and the distribution of the 
moments is Gaussian: 



(X)=Q,(X 2 ) = l,(X 3 )=Q, (X 4 ) = 3, . 



(37) 



The notion that it is impossible to achieve directed mo- 
tion (or equivalently, to extract work) from a system in 
thermal equilibrium is confirmed. At least two reservoirs 
at different temperature are necessary to break detailed 
balance and make possible the rectification of thermal 
fluctuations. 



C. Asymmetry 

The geometry of the motor units enters into the ex- 
pression of the average speed via de shape factors <ji and 
(j 3 . In table [i] we have reproduced these quantities for 
the objects depicted in Fig. [3] As is clear from symme- 
try arguments, the appearance of systematic motion in 
one direction requires, apart from non-equilibrium con- 
ditions, also the breaking of the spatial symmetry in the 



system. One easily verifies from Eq. ( 30 1 that 
Vn,i = 0, n odd, 



(38) 



when the surface of a motor element possesses reflection 
symmetry along the z-axis, the direction of motion. Con- 
sistent with this symmetry observation we find that the 
drift speed of the motor in steady state is indeed zero 
at lowest order in the perturbation when 03^ = 0, cf. 
Eq. (351. It is however interesting to note that reflection 



symmetry is a sufficient but not a necessary condition for 
a 3,i — 5 an d hence for obtaining zero sustained motion 
(at least in this order of the approximation). Consider 
for example a spherical cone, see Fig. [5] For a specific 
opening angle of « 55° 03^ becomes zero, even though 
there is no reflection symmetry. A similar discussion can 
be applied to higher orders corrections in the e-expansion 
(featuring the appearance of the higher shape factors 175, 
CT7, and so on) indicating that there are special shapes 
which will have a very low average speed even though 
there are no immediate symmetry reasons to expect so. 
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FIG. 5: The geometric factors ai/S and oz/S as a function 
of the opening angle a for a spherical cone, with S the sur- 
face area. While 02 appears in the expression for the friction 
coefficient, (73 relates to the drift speed of the motor. When 
the geometry of the surface exhibits symmetry along the z- 
axis (the direction of movement), 03 is zero, and the motor 
shows no directed motion. The spherical cone is an example 
of a class of shapes for which 03 can become zero, and hence 
the drift speed (at least to the first approximation), with- 
out however showing reflection symmetry. For this particular 
case, 03 = for a m 55°. 
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FIG. 6: The drift speed as a function of the temperature 
difference AT = Ti — T2 between the reservoirs for two con- 
figurations with identical motor units: parallel (solid curves 
and lower inset) and antiparallel (dashed curves and upper in- 
set). The units are cone-shaped (opening angle 30°) and the 
motor mass is 1000 kDa. Both reservoirs are filled with argon 
gas of same density and T2 is fixed at 299 K. For small AT 
(see insets), the dependence is parabolic for the parallel setup 
and linear for the antiparallel setup. The drift speed is of 
order nm/s (parallel) and /im/s (antiparallel) for AT m IK. 



D. Temperature gradient 

For simplicity we limit the discussion of the drift ve- 
locity, in particular in relation to the applied tempera- 
ture gradient, to the case of two reservoirs. Furthermore, 
the impact of the geometry is most clearly demonstrated 
when the motor units have an identical shape in both 
reservoirs. There are then two possibilities, namely either 
the units have the same orientation (parallel) , or they are 
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pointing in opposite direction (antiparallel) , see Fig.[6]for 
a schematic representation. For the first scenario, with 
C2.1 = °2,2 = <?i and a 3i i = a 3i2 = cr 3 , 



(V)e = 



Kk B P\Pi(T\ 



Eq. (|35|) yields: 



8M 



(piVTi + P2\ r T 2 )' £ 



a-. 



(39) 

For the second scenario, careful consideration of the sign 
of cos# in Eq. (30) leads us to write cr 2 ,i = 0-2,2 = 02 as 
in the first scenario, but now when cr 31 = a 3 , it follows 
that (T3.2 = — 03. The expression for the drift speed thus 
becomes: 



(V)o = 



nk B pxp2@\ 
8M ( m 



-r 2 )(V7Y + v 



r 2 ) 0-3 



(T 2 



One striking feature of these results is that the drift ve- 
locity is scale-invariant. This is a general property: both 
<j 2 and 03 scale linearly with the total surface S of the 
motor units. As they appear in the denominator and 
the nominator respectively in Eqs. (39 40 1, the scale de- 
pendence cancels out. This becomes even more appar- 
ent for the particular cases presented in Table [I] where 
<j 2 /S and cr 3 /S are expressed in topological terms. Note 
that the scale invariance of the drift speed is only valid 
with respect to the scale of the entire motor. The rela- 
tive proportions of separate motor units do matter. Note 
also that scale invariance applies when disregarding the 
dependence on the mass M. For comparison with a phys- 
ically realistic situation, we assume a constant density of 
the motor, so that the drift velocity will decrease with in- 
creasing size of an object, through its 1/Af-dependence. 

To investigate the departure from the equilibrium 
state, we consider small deviations of T\ and T 2 about 
the average T, 



T X = T 



AT 

~2~' 



T 2 = T- 



AT 



(41) 



so that for AT/T < 1 the drift speeds tend to: 



{V) t 



(V)o 



1 

16 



P1P2 C3 
M (pi + p 2 ) 2 a 2 



k B T (AT\ 
M \ T ) 



P1P2 cr 3 
M (pi + p 2 ) 2 a 2 



k B T ( AT\ 
M \ T ) 



(42) 
(43) 



We note that the respective orientation of the motor 
elements, parallel or antiparallel, plays a crucial role. In 
the first case the sustained average displacement is al- 
ways in the same direction (that of — 03), indifferent of 
the sign of AT. This was to be expected because there is 
an additional symmetry in the system: the interchange 
of the temperature reservoirs has no effect. From the 
point of view of irreversible thermodynamics, this situ- 
ation is special since there is no linear relation between 
the thermodynamic force (the temperature gradient) and 
flux (the resulting speed of the motor). For the second 



TABLE II: Values for the drift speed of the motor, as pre- 
dicted by theory. Realizations of the motor at different length 
scales from micrometers to nanometers are presented. The 
motor units are modeled as silica (SiCh) cone-shaped objects, 
located in two reservoirs containing argon gas at temperatures 
299.0 K and 299.1 K. For constant motor mass, the speed in- 
creases as the opening angle a of the cone decreases. 



a = 30° 



a = 5° 



a = 30° 



a = 5° 



Motor 

mass 

(kDa) 


Cylinder base 
(nm) 


Drift speed 
(pm/s) 




10 iu 


2600 


1400 


5.2xl0" s 


9.5x10" 


-8 


10 7 


260 


140 


5.2xl0~ 5 


9.5x10" 


-5 


10 6 


120 


55 


5.2xl0~ 4 


9.5x10" 


-4 


10 s 


63 


29 


5.2xl0" 3 


9.5x10" 


-3 


10 4 


26 


14 


0.052 


0.095 




1000 


12 


6.3 


0.52 


0.95 




100 


5.5 


2.9 


5.2 


9.5 




10 


2.6 


1.4 


52 


95 





case of antiparallel alignment, one observes the usual sit- 
uation of linear response between thermodynamic force 
and flux [21 E]: equilibrium is a point of flux reversal, the 
direction of net motion reversing with AT-inversion. 

As an illustration, we reproduce, in table [IT] explicit 
values for the drift speed in the case of a single asym- 
metric unit, namely a cylindrically symmetric cone, po- 
sitioned in antiparallel alignment in the two reservoirs 
under physically realistic conditions. The degree of asym- 
metry is described by a single parameter, the opening 
angle a. The geometric factors a 2 and a 3 are known an- 
alytically (see Table |l| and we find the following simple 
expression for the drift speed: 



(V) = 




Ttk B pip 2 (Ti - r 2 )(v7T 
8M (piVT~i + p 2 VT 2 ~) 



2 (sin a- 



!)• 

The drift speed will become zero for a = 90°, namely 
when the cone loses its asymmetry and reduces to a flat 
disk. A natural question is whether there is an optimal 
opening angle a that maximizes the drift speed. For a 

fixed cross section, one finds a = scc _1 [-y/ (1 + v / 5)/2] ~ 
38°. If, on the other hand, one assumes that the mass 
is kept constant, a maximal speed is reached for an in- 
finitely sharp cone. 

Note finally the very strong size-dependence: objects 
of 20 nm cover their length 5 times per second, for 5 nm 
size objects this becomes 1200 times per second. 



VI. CONCLUSION 

We have calculated, on the basis of an exact micro- 
scopic theory, the properties of a thermal Brownian mo- 
tor in a three-dimensional setup. When detailed balance 
is broken by the application of a temperature gradient, a 
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systematic net speed appears, as given in Eq. (35 1. As 



an example, for a motor consisting of cone-shaped silica 
units of size 20 nm, one obtains a drift speed of about 0.1 
jim/s when subject to 0.1 K temperature difference in a 
gaseous environment. It remains to be seen whether the 



predictions of our theoretical analysis (involving various 
simplifications such as molecular chaos, elastic and nor- 
mal interactions between gas particles and motor, expan- 
sion in mass ratio) provide a realistic estimate, especially 
for motors operating in a viscous environment. 
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